Electronic structure and magnetic properties of pyroxenes (Li,Na)TM(Si,Ge)20 

Low-dimensional magnets with 90° bonds 
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The results of the LSDA+U calculations for pyroxenes with diverse magnetic properties 
(Li,Na)TM(Si,Ge)206, where TM is the transition metal ion (Ti,V,Cr,Mn,Fe), are presented. We 
show that the anisotropic orbital ordering results in the spin-gap formation in NaTiSi2 06. The 
detailed analysis of different contributions to the intrachain exchange interactions for pyroxenes is 
performed both analytically using perturbation theory and basing on the results of the band struc- 
ture calculations. The antiferromagnetic t^g — t^g exchange is found to decrease gradually in going 
from Ti to Fe. It turns out to be nearly compensated by ferromagnetic interaction between half-filled 
t2g and empty Cg orbitals in Cr-based pyroxenes. The fine-tuning of the interaction parameters by 
the crystal structure results in the ferromagnetism for NaCrGe2 06. Further increase of the total 
number of electrons and occupation of eg sub-shell makes the t2g — eg contribution and total ex- 
change interaction antiferromagnetic for Mn- and Fe-based pyroxenes. Strong oxygen polarization 
was found in Fe-based pyroxenes. It is shown that this effect leads to a considerable reduction of 
antiferromagnetic intrachain exchange. The obtained results may serve as a basis for the analysis 
of diverse magnetic properties of pyroxenes, including those with recently discovered multiferroic 
behavior. 

PACS numbers: 72.80.Ga, 71.20.-b, 74.25. Ha 



I. INTRODUCTION 

Pyroxenes with the general formula AMX2O6 are an 
extremely rich class of compounds. Here A may be alkali 
Na, Li or alkaline-earths Ca, Sr,... elements, M- different 
metals with the valence 3-|-, X is most typically Si^"'', but 
also Ge^+ and even could be . These systems are best 
known in geology and mineralogy: they are one of the 
main rock- forming minerals in the Earth's crust, and are 
also found on the Moon, planets and in meteorites. One 
of the pyroxenes, NaCrSi206, may be found in natural 
form only in meteorites, which is reflected in its name: 
mineral kosmochlore. Another one, NaAlSi2 06, is the 
famous Chinese jade. For the condensed matter physics 
the pyroxenes with M a transition metal (TM), are of 
particular interest because of their rich and interesting 
magnetic properties. 

Structurally magnetic sub-system in pyroxenes is 
quasi-one-dimensional: TMOg octahedra form edge- 
sharing zig-zag chains connected by the chains of XO4 
tetrahedra (see Fig. [1]). Two aspects are here especially 
important. First of all the edge-sharing character with 
TM -O- TM angle close, but not exactly equals to 90°, 
which results in the competition between different contri- 
butions to the super-exchange (see Sec. HI]). Secondly a 
"shifted" packing of neighboring TM-chains leads to the 
situation, when every TM ion is connected (via own O 
and XO4 tetrahedra) with two TM ions in neighboring 



chain. Thus, the general topology of exchange interac- 
tion is triangular- like, see Fig. HJa), and the magnetic 
system may be frustrated. 



Both these factors may lead to nontrivial magnetic 
properties of pyroxenes, which are indeed observed exper- 
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FIG. 1: (color online) Crystal structure of pyroxenes. Two 
chains made of TMOg octahedra and its connections via XO4 
(where X is usually Si or Ge) tetrahedra are shown. Purple 
balls represent Li(Na) ions. 
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FIG. 2: (color online), (a) Triangular- like topology of the 
TM network in pyroxenes. Solid lines are intrachain links, 
dashed and double-dashed lines are the interchain connection 
via one or two XO4 tetrahedra; J, Ji and J2 are correspond- 
ing exchange parameters, (b) Sketch of the one of TMOg 
chains. Two different types oxygens are shown. Those de- 
picted as solid circles belongs to two different TMOe octahe- 
dra, open circle - only to one. The later oxygens are called 
"side-oxygens" in the text. The natural choice of the local 
coordinate system (LCS), when 2:— axis is directed along the 
single 180° 0-TM-O bond connecting two neighboring TMOe 
octahedra, is shown. 



imentally, although for most of them the detailed studies 
arc still absent. There are among them antiferromag- 
netic (e.g. (Li,Na)V(Si,Ge)206 H, (Li,Na)FeSi206 M) 
and even ferromagnetic (e.g. NaCrGe206 0) insulators. 
Some of the pyroxenes have the ground state with the 
spin gap: CaCuGeaOg 0, NaTiSiaOg Si, although 
the nature of such a state is apparently different: this 
seems to be due to formation of interchain dimers in 
CaCuGeaOe 0, and due to the intrachain dimerization 
in NaTiSi20^ stabiHzed by the special type of orbital 
order 

Finally a new twist in this story is the recent dis- 
covery that the pyroxenes apparently form a new class 
of multiferroics [Ij. It was found that at least three of 
them, NaFeSi206, LiFeSi2 06 and LiCrSi206, develop an 
electric polarization in a magnetically-ordered state, be- 
low - 6 K in NaFeSi2 06, ~ 18 K in LiFeSi2 06 and 
11 K in LiFeSi2 06- Similar to other multiferroics 
with magnetically-driven ferroelectricity (see e.g. re- 
views [HjSI), electric polarization in them can be strongly 
modified by a magnetic field. Note also that the most 
pyroxenes are insulating (often beautifully-colored crys- 
tals), which makes them especially good for studying 
eventual ferroelectricity. 



To understand magnetic properties of this large and 
important class of compounds, and especially to have a 
better understanding of the possible mechanism of mul- 
tiferroic behavior in them, it is desirable to know the 
details of the exchange coupling. The knowledge of or- 
bital occupancy is also very important, both in connec- 
tion with their structure and for determining the details 



of exchange interactions. To reach this goal, and also to 
provide better understanding of an electronic structure 
of pyroxenes in general, we undertook a detailed theo- 
retical investigation of these problems, concentrating on 
systems with A=Na, Li; X=Si,Ge and different trivalent 
TM ions. We carried out detailed ab-initio calculations, 
and in particular obtained the values of the exchange 
constants. Besides these values themselves, which can 
be useful for interpreting the properties of real materials, 
we paid main attention to the study of the systematics of 
the exchange in different compounds of this large class, 
and to particular mechanisms, or partial contributions 
of different exchange passes to the total exchange. Gen- 
eral conclusions reached are important for the analysis of 
the multiferroic behavior in pyroxenes, they may be use- 
ful also for understanding of the properties of pyroxenes 
with alkaline-earths A ions such as e.g. CaMnSi206 and 
may serve as a very good illustration of different tenden- 
cies in exchange couplings in systems with complicated 
geometry in general. 



II. DIFFERENT CONTRIBUTIONS TO 
EXCHANGE 

Before presenting the results of our ab-initio calcula- 
tions, let us first discuss different mechanisms of the in- 
trachain exchange interaction. The explicit expressions 
for it may be calculated using the perturbation theory in 
t/Udd or t/ A, where t is an effective electron hopping (ei- 
ther direct d~d hopping tdd or hopping between d— states 
of TM and 2p— states of oxygen tpd), Udd is an on-site 
Coulomb interaction, and A is a charge-transfer energy 
(energy of the promotion of an electron from 2p— states 
of oxygen to Sd— states of TM). 



The main motive of the intrachain packing of TM 
ions in pyroxenes is the edge-sharing packing of TMOg 
octahedra (see Fig. [T]), with 90° TM-O-TM bonds for 
(half)filled t2g — t2g and <2g — Cg orbitals. This implies 
that one of the important contributions will be a direct 
exchange between partially- filled t2g orbitals, the lobes of 
which look towards neighboring TM ions. This contribu- 
tion to the total exchange will be AFM and proportional 
to 

J t-2g-t-2g /-. S 

•-'t^g-t^g ~ Udd ' 

where tt2g-t2g is the t2g ~ t2g hopping parameter due 
to direct overlap of d— wave functions. Since consecu- 
tive TM-O-TM planes in zig-zag TM-chains are rotated 
by 90° (see Fig. [2b), specific orbital pattern will be cru- 
cial for this contribution. This, in particular, leads to 
strongly alternating exchange in NaTiSi206, determin- 
ing the appearance of the spin gap in it (see Sec. lIVp . 
On the other hand, corresponding orbital occupation in 
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FIG. 3: (color online). Typical orbital ordering between two 
occupied TM d— orbitals via the same O— 2p orbital, which 
results to AFM interaction for the case of 90° TM -O- TM 
bond. 




FIG. 4: (color online) . Typical orbital ordering between occu- 
pied (filled) and unoccupied (unfilled) d— orbitals of TM via 
the same 0—2p orbital, which results to FM interaction for 
the case of 90° TM -O- TM bond. 



V pyroxenes does not lead to an alternating exchange 
(Sec. IYJ. And of course for occupation, as for Cr^+ 
or Fc'^+; all the pairs between neighboring TM ions be- 
longing to the same chain would have equal contributions 
to this part of the exchange. 

Due to decrease of the radius of d— orbitals and with 
it of t2g — t2g hopping, this direct exchange is expected 
to strongly diminish in going from Ti to Fe pyroxenes, 
see below. Note that we use the terminology "direct ex- 
change" for the exchange caused by the direct d — d over- 
lap and hopping, and reserve the term "super-exchange" 
for the exchange mediated by oxygens [l^ . 

Besides direct d — d hopping, p — d hybridization leads 
to an important contribution to the total exchange due to 
TM-O-TM hopping. Different possible situations leading 
to a super-exchange via oxygens are illustrated in Figs.[3]- 
[5j In these figures (half) occupied orbitals of TM are 
shown as filled (blue or hatched red) , and empty orbitals - 
as empty (white) ones. ^From these figures it is clear that 
the corresponding contributions can be both antiferro- 
and ferromagnetic. 

The simplest case is that of Fig. EJa). Here half- filled 
orbitals at neighboring TM's overlap with the same oxy- 
gen p— orbital. According to the usual rules [l^ this 
antiferro-orbital ordering gives a substantial AFM ex- 
change: 

jSE ''pd.TT ( ^ J_\ 

I2ATIV V^i' ' 

here tp^Tr is the tt— hopping between tig and 2p— orbitals. 
It is important that the magnitude of tig — O — tig {eg — 
O — Cg) exchange presented in Eq. ([2|) does not strongly 
depend on the TM-O-TM angle a (only via dependence 



of tp^ hopping on a). First and second contributions 
are due to the correlation and delocalization effects in 
terminology of Ref . [isl. 

Similarly, in the case of Fig. Elb) half-filled tig and eg 
orbitals overlap via the same oxygen orbital, which gives 
even stronger AFM exchange 

where tpda is the hopping between eg orbital of TM and 
2p— orbital of oxygen directed towards it. tpda is approx- 
imately two times larger than tpd-w [T3 | . 



If the hopping is only possible from the occupied to an 
empty d— orbital, or if occupied orbitals of TM do not 
overlap with the same oxygen orbital, the resulting con- 
tribution to the exchange turns out to be ferromagnetic. 

First of all let us consider the situation when one of 
d— orbitals is empty and it overlaps with the same oxygen 

orbitals, as the occupied d— orbital, see Fig. Hl^a.b). 
Then electrons can virtually hop to this empty orbital, 
but the Hund's rule coupling J^^^ tends to orient spins 
at an ion parallel. This will result in a FM exchange: 

jSE ''pdrn'-pdm' ( -J h , H \ 1 a\ 

•^t2,-t2,/e, ~ A2 V (2A + C/pp)2 ^ / ■ ^' 

Note that J^*^ here is the Hund's energy gain when we 
put an extra electron with parallel spin to the original 
state of a TM ion, i.e. depending on the total spin of the 
latter J™ may be equal e.g. to twice the usual atomic 
value of Jh '^0.8 eV for 3d— series when the spin of the 
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FIG. 5: (color online). Typical orbital ordering between two 
occupied d— orbitals of TM via different O— 2p orbitals, which 
results in a FM interaction for the case of 90° TM -O- TM 
bond. 



TM ion is S = 1 (two electrons with parallel spins as 
in V3+), or 3Jh, as for Cr^+ {(f, S = 3/2). This may 
enhance the ferromagnetic contribution in certain cases 
(in particular for Cr systems). 

The last contribution to the total d — d exchange inter- 
action is coming from the overlap of occupied rf— orbitals 
with different (orthogonal) 2p— orbitals. Typical orbital 
pattern for this situation is shown in Fig. [5l but one 
may draw similar pictures also for the case of eg/sg or 
^g/hg orbitals. This exchange will be FM and again be- 
cause of the Hund's rule, but in this case acting in oxygen 
2p— shell. With such an overlap the consecutive hopping 
of the electrons from one d— ion to another one through 
oxygen p— shell (delocalization exchange of Ref. [Tsl ) turns 
out to be forbidden. Two oxygen's electrons are hopping 
to different TM ions instead. Because of the Hund's cou- 
pling on oxygen (Jfj) the FM ordering becomes more 
favorable: 



jSE 



-tigje 



f2 ^2 TP 

"A2(2A + C/pp)2 



(5) 



Note that tpdm in the expressions ([4]) and ([5]) may be tpd-K 
or tpda, depending on the particular case. For instance, 
for the case depicted in Fig. IHb) one them is iptj^r and 
another tpda- 

In order to describe magnetic interactions in a particu- 
lar pyroxene one has to include in general all possible con- 
tributions discussed above, via both common oxygens, 
and take into account several general trends. 



First of all, the radius of d— states (r^) and with them 
the direct d — d (tdd) hoppings decrease going from the 
light to heavy TM (from Ti, V to Fe). According to the 
famous Harrison scaling , 



tdd 



' d 



(6) 



where D is the distance between TM. Taking into account 
typical values for presented in Ref. [13 and assuming 



only this effect may provide suppression of the direct 
exchange in Fe-based pyroxenes comparing with those 
based on Ti by a factor of 6. Secondly, the distance be- 
tween TM ions D by itself strongly affects the exchange 
coupling, as seen from Eq. ([6]). Third, there exist the 
strong increase of Udd and decrease of the charge-transfer 
energy A for the late TM's. A enters the expressions 
Eq. (HI)- ([5]) in the denominators either as A^ or even as 
A'^ and plays very important role. For Ti'^+ A equals to 
6 eV, whereas for Fe'^+ it is 2.5 eV only [l^. This factor 
leads for the late 3d— elements to an increase of the role 
of super-exchange via oxygen, which competes with the 
direct d — d exchange dominating in early TM's. 

The qualitative discussion presented in this section will 
serve as a framework for interpreting the results of the 
following sections. As we will see, all these factors are 
quite important and indeed determine the tendencies of 
exchange constants in pyroxenes and many other systems 
with 90° TM -O- TM bonds. 



III. COMPUTATIONAL DETAILS 

Crystallographic data used in the calculations 
were taken from the follo wing papers: NaTiSi206 
(T==100 K, space group P-1) [Tg, LiVSizOg (T=293 K, 
space group C2/c) [H, NaVSizOg (T=296 K, space 
group C2/c) [13, LiCrSizOg (T=100 K, space group 
P2i/c) Jii, NaCrSizOg (T=300 K, space 
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23, NaFeSizOe (T=14 K, space group C2/c) [2 



that D is the same, we found that 



2-5i£f. Thus, 



C2/c) 
C2/c) 
C2/c) 

LiFeSizOg (T=100 K, space group P2i/c) 

There is a controversy in the literature about the crys- 
tal structure of NaVGe206 and LiVGe206. Although 
all the studies were performed for the room tempera- 
ture, the difference in lattice parameters exceeds 0.3 A 
for NaVGe206 [3, [ISl- There is also uncertainty about 
the space group for this material [l^, [IH] . In contrast to 
NaVGe206, there is no striking disagreement about the 
crystal structure of LiVGe206 in Ref .13. [2^ However, we 
were not able to reproduce distances presented in Ref. [23 
from their crystal structure. Therefore, the study of elec- 
tronic and magnetic properties of both NaVGe206 and 
LiVGe206 was postponed for the future, when all the 
problems about the crystal structure of these compounds 
would be resolved. 

The primary calculations were performed within the 
framework of the linear muffin-tin orbitals method [27t . 
The values of on-cite Coulomb interaction ([/) and 
Hund's rule coupling {Jh) parameters were taken as fol- 
lowing: Uti = 3.3 eV, Uv = 3.5 eV, Ucr = 3.7 eV, 
UMn,Fe = 4.5 cV; J H = 0.8 eV for Ti, V and Cr, 
Jh == 0.9 eV for Mn and = 1 eV for Fe. The values 
of U and Jh for Ti were calculated for the same Ti radii 
in Ref. [1^. For other TM those parameters were taken 
as an average in wide variety of values of U and Jh for 
TM'^^ ions presented in the literature. For the material 
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with small values of exchange parameters (cremates) we 
checked that the small variation of U {±0.5eV) doesn't 
change magnetic ground state. 

Instead of LDA+U approximation used in Ref. we 
utilized LSDA+U approach in the present paper. The 
difference between these two methods is in the parts of 
the Hamiltonian allowed to be spin polarized. In the 
LDA+U method only d-shell of the transition metal ions 
may be spin-polarized; for all the other ions different 
spin-states are averaged out. This allows more easy ac- 
count of the double counting in LDA-I-U formalism, but 
does not treat correctly the effects of possible spin polar- 
ization of oxygens, which significantly restricts the ap- 
plicability of the method to more complex systems. We 
found the noticeable spin polarization of oxygens for the 
late TM pyroxenes (Mn, Fe, see SedVlII in the LSDA+U 
calculations and used this approach for the whole se- 
ries. For the Ti and V systems the effects of oxygen po- 
larization is insignificant and these two approximations 
give essentially the same results. For Fe-based pyroxenes 
LDA+U overestimated exchange constants by as much 
as ~ 25 - 30%. 

The second difference from the Ref. is the more care- 
ful choice of the MT-radii. We found that for not uniform 
atom distribution in dimerized NaTiSi2 06 the LMTO 
calculation is very sensitive to the radii of Na MT-sphere. 
As a result of too large RNa — 3.5 a.u. used in Ref. 
Na ion takes the electrons in excess, which leads to the 
significant overestimate of the exchange constant. In the 
present paper we used Rnu = 3.3 a.u. for all pyroxenes. 
This choice of RNa results in a good agreement with the 
band structure and exchange constants obtained in the 
pseudo-potential calculations, as it will be discussed in 
Sec. IIVI Note, however, that the main conclusions about 
the importance of the strong Coulomb correlations made 
in Ref. are valid, and do not depend on the choice 
of RNa- In addition, it is interesting enough that only 
the results for NaTiSi206 significantly depend on the 
RNa- Other pyroxenes do not show such a strong de- 
pendence. For instance, by going from Rnu = 3.5 a.u. 
to RNa = 3.3 a.u. for NaCrGe206, the exchange con- 
stants change by 1.1 K, and for NaCrSi206 - even less, 
0.2 K. 



The radii of other elements were set as following: Rq = 
1.65 - 1.75 a.u., Rsz,Ge = 1-9 a.u., Rlz = 2.9 a.u. The 
radii of transition metals were chosen to be Rtm ~ 2.4 — 
2.67 a.u. 

The Ti,V,Cr,Mn,Fe(4s,4p,3d), 0(2.s,2p,3d), 
Si(3s,3p,3d), Ge(4s,4p,4d), Li(2s,2p,3d) and 
Na(3s,3p,3c?) orbitals were included to the orbital 
basis set in our calculations. The Brillouin-zone (BZ) 
integration in the course of the self-consistency it- 
erations was performed over a mesh of 64 k-points 
in the irreducible part of the BZ. For those systems 
for which exchange parameters are particularly small, 
the calculation was checked by the finer mesh of 216 
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FIG. 6; (color online). Total and projected DOS for 
NaTiSi2 06 obtained in LSDA+U calculation. Total DOS is 
the sum for both spins. Up (down) panel corresponds to spin 
majority (minority). Fermi level is zero energy. 

k-points. 

In our notations, Heisenberg Hamiltonian is written in 
the following form: 

H = Y,J,,Sj,, (7) 

ij 

where summation runs twice over every pair J were 
computed in the framework of Lichtenstein's exchange 
interaction parameter (LEIP) calculation procedure [28| . 
According to the LEIP formalism, exchange constants J 
can be calculated as the second derivatives of the energy 
variation at small spin rotation. For s = 1/2 

^«J = V /* ^ A™"'G™'f" Af ""'g;^7™, (8) 

-oo TJ™',,, 
m m 

where m, m' , m", m'" are the magnetic quantum num- 
bers, i,j - the lattice indexes, Af""' = H^f-H^f is 
the on-site potential correction due to different spins and 
the Green's function is calculated in the following way: 

„mn m'n * ( l.\ 

k,n " 

Here a is the spin index, c™" is a component of the n-th 
eigenstate, and E^ is the corresponding eigenvalue. Since 
transition metal ions are usually magnetic, the Green 
functions G-"™ (e) and potential corrections A™™ are 
computed only for them. 

The great advantage of the LEIP procedure is the pos- 
sibility to calculate not only total J's between corre- 
sponding sublattices of transition metal ions, but, since 
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FIG. 7: (color online). The orbital ordering obtained in 
LSDA+U calculations for NaTiSi2 06 is shown. 



it is linear with respect to the magnetic quantum num- 
bers, one can also compute partial contributions, which 
come from different orbitals. 

These orbitals were defined in the local coordination 
system (LCS), where axes are orthogonal and directed 
as much as possible to the neighboring oxygens. The 
natural choice of the LCS is when the local z— axis is 
directed along the only 180° O — TM — O connecting 
two neighboring TMOg octahedra, see Fig. [Hb). This 
definition is unique, and it is used for all the pyroxenes 
studied in the present paper. 

Since some of the exchange constants obtained within 
the LMTO method using LEIP were found to be sensi- 
tive to the choice of MT-radii, in the present paper we 
recheckcd the LMTO results by calculating the exchange 
constants from the total energies in the pseudo-potential 
PW-SCF code [1^. The formalism of pseudo-potential 
method does not require the introduction of any atomic 
spheres. We used ultrasoft pseudo-potentials with non- 
linear core correction taken from www.pwscf.org with dif- 
ferent forms of exchange-correlation. The plane- wave and 
kinetic energy cut-off's were chosen to be 40 Ry and 200 
Ry, respectively. 216 A:— points were used in the course 
of self-consistency. 

Hopping integrals were obtained using Wannicr func- 
tion projection in LMTO method described else- 
where [l^, [S^l . Small subspace defined by three t2g or- 
bitals (defined in LCS) per V atom was used in this pro- 
cedure. Since for all of the considered pyroxenes t2g- 
bands are well isolated, resulting projected and initial 
LDA bands completely coincide. 



IV. Ti-BASED PYROXENE: NaTiSiaOs {d}) 

As it was discussed in details in Sec. U pyroxenes 
are made of isolated zig-zag chains of TMOg octahe- 
dra. Combination of low dimensionality together with 
the small value of the spin moment in NaTiSi206 results 
in an opening of the spin-gap [1, [13 . Here we will show 
that it becomes possible due to a particular orbital or- 
dering. 



Two different models were suggested in order to ex- 
plain the physical origin of the spin-gap formation. The 
first one (dimer model) was initially proposed in exper- 
imental papers [1, [3l| and supported afterward by the 
LDA+U calculations Q. It takes into account quite 
strong dimerization of Ti-Ti chain found in very de- 
tailed structural study performed by Redhammer and 
co-authors [l^. According to this model the ground 
state of NaTiSi206 consists of spin singlet Ti-Ti dimers, 
and the spin-gap observed in the experiment is the en- 
ergy gap between spin singlet and triplet states of these 
dimers. 

In another model (Haldane scenario) [H, [s^l Ti-Ti 
dimers were proposed to be in a spin triplet state and 
such S* = 1 Ti-Ti pairs form the AFM chain, which is 
known as Haldane chain [s^ with the spin-gap l^sG = 
0.41J, where J is the exchange constant. 

Thus, the key point in these two scenarios is the de- 
generacy of the magnetic ground state of Ti-Ti dimer: in 
the dimers model it is a singlet, but in the Haldane one 
- a triplet. 

It is interesting to note that depending on the situa- 
tion the pair of metal ions may have quite different de- 
generacy of the magnetic ground state. In the case of 
isolated cluster of metal atoms in vacuum, the ions in 
the pair are very tightly bound and must be considered 
as one quantum-mechanical object. The bond length of 
such isolated "metal-metal molecule" respectively is quite 
small: 1.9 A for Ti-Ti pair [3^1. The spin configuration in 
this case will be governed by the Hund's coupling, which 
favors the state with larger multiplicity. As a result the 
ground state of isolated metal-metal dimers usually is not 
the singlet one [s^. 

However the situation in insulating transition metal 
oxides is qualitatively different. Typical metal-metal dis- 
tance varies here depending on the type of the crystal 
structure, but usually it is much larger than in isolated 
dimers: Ti-Ti distance is 3.1 A in the edge-sharing oc- 
tahedral structure of NaTiSi2 06 [Q. As a result every 
ion has its own ground state configuration, and the in- 
fiuence of the other ions can be incorporated via per- 
turbation theory [l^. The resulting electronic Hamilto- 
nian can be reduced to a Heisenberg model (36| . which 
is known to have spin-singlet ground state for the dimer. 
Our band-structure calculations in LSDA+U approxima- 
tion give the same result: rather strong antiferromagnetic 
exchange in short Ti-Ti dimers Jintra — 396 K, and the 
weak exchange between dimers Jinter = 5 K. 

The electronic structure of NaTiSi2 06, obtained in the 
LSDA+U approach is presented in Fig. [6] Both the 
top of the valence band and bottom of conduction one 
are formed by Ti-Sd states, divided by the band gap of 
1.77 eV. Singly occupied Ti-3c? orbital has zy— character. 
Corresponding ferro-orbital ordering is shown in Fig. [T) 
In a short Ti-Ti pair it results in a strong direct ex- 
change interaction, which is AFM in agreement with 
Goodenough-Kanamori- Anderson (GKA) rules [l^, . 
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On the other hand, for the long Ti-Ti pairs the orbitals 
are directed in such a way that the overlap between them 
is almost zero, which leads to a rather small inter-dimer 
exchange. 

Note that the exchange constantspresented above are 
different from those given in Ref. |9|. This is due to a 
more careful choice of Na MT-radii, as it was discussed 
in Sec. IIIII In the present paper we additionally checked 
the LMTO results by calculating total energies using 
pseudo-potential PW-SCF code [2^. In order to have 
reliable result the calculations with the use of two differ- 
ent pseudo-potentials based on Perdew-Wang-91 (PW91) 
and Perdew-Bruke-Ernzerhof(PBE) exchange potentials 
were utilized. The total energy difference between FM 
and AFM solutions were found to be 398 K for PW91 
and 395 K for PBE, which arc very closed to the present 
LMTO results. 

One should be careful comparing exchange constants 
obtained in the band structure calculation with the value 
of the spin gap. This is due to strong quantum effects. 
What is indeed calculated by any mean-field method are 
the total energies for AFM and FM configurations (or the 
change of these energies with respect to the small rotation 
angle of spins, as in LEIP, which is essentially the same). 
For the case of isolated dimer s = 1/2 and Heisenberg 
model defined by Eq. ([7]) the total energy difference be- 
tween FM and AFM configuration equals to J. The spin 
gap AsG is the energy difference between quantum sin- 
glet and triplet states, which are not the same as coUinear 
AFM and FM states. For the 5 = 1/2 the spin-gap 
AsG = 2 J. The spin-gap obtained for NaTiSi206 using 
different experimental methods is about 600-700 K [STj . 
This agrees with Asg ~ 790 K calculated in the pre- 
sented paper. 

V. V-BASED PYROXENES {(f) 

Ionic configuration of V in the pyroxenes is (P imply- 
ing S = 1. Because of a larger spin moment it is much 
harder to make dimerized chains in V-based pyroxenes. 
The spin gap may still appear in the system, if it could 
be considered as a set of isolated one-dimensional AFM 
chains of 5" = 1 (Haldane chains [l^l). However, exper- 
imentally all the V-pyroxenes show long-range magnetic 
ordering at low temperatures (3]. Three different rea- 
sons for it were discussed in the literature: biquadratic 
exchange [H, [s^ , interchain exchange [iO, |4l| and next- 
nearest- neighbor coupling [ssj . 

Common structural feature for all of the pyroxenes ex- 
cept Mn-based ones (they are Jahn- Teller active) is the 
presence of two short TM-0 bonds. These are the bonds 
with "side-oxygens", see Fig. ??(b). 

Short V-0 bonds form nearly 90° angle. Such a dis- 
tortion of TMOg octahedra may be represented as two 
independent compressions in x and y directions. As a re- 
sult zx— and zy— orbitals will be almost degenerate and 
have lower energy than xy. 



TABLE I: Results of the LSDA-I-U calculations for different 
pyroxenes. The values of the band gap and spin moments were 
taken from AFM (intrachain) solution intrachain calculations 
for all the compounds, except NaCrGe2 06. The exchange 
parameters were obtained in LEIP formalism 



Ionic conf. Band Spin Exchange 
of TM gap (eV) mom. (/-ts) intrachain (K) 

NaTiSi2 0c 1.77 0.92 396 

LiVSi206 (P 1.92 1.87 50.3 

NaVSi206 1.80 1.85 15.0 

LiCrSi206 d'^ 3.57 2.78 3.4 

NaCrSi206 d^ 3.03 2.83 -0.8 

NaCrGe2Q6 d^ 2.88 2.84 -5.2 

NaMnSi2 06 rf^ 1.76 3.69 3.4 

LiFeSi206 d^ 1.63 4.02 15.8" 

NaFeSi206 d^ 2.31 4.24 15.9° 



"Note that because of the large O— 2p polarization the exchange 
parameters obtained in LEIP formalism are significantly overesti- 
mated for Fe-based pyroxenes; see Sec. IVIIll 

In our LSDA-I-U calculations two 3d electrons of V 
follow the crystal-field splitting and occupy zx— and 
zy— orbitals, i.e. these systems do not have orbital de- 
generacy. Each orbital provides strong AFM direct ex- 
change with one of the neighboring V. This leads, as we 
will show below, to the uniform exchange interaction in 
each chain, with the direct exchange defined in Eq. ([1]). 
The super-exchange via oxygens would be also homoge- 
neous and (weaker) FM. 

In our calculations all the (Li,Na)VSi206 pyroxenes 
were found to be insulators with the band gap 1.8 — 
1.9 eV. Spin moments are slightly reduced from ionic 
value of 2hb because of the hybridization. Interestingly, 
the calculated intrachain exchange constants significantly 
decrease going from Li i— > Na (see Tab. |l]). In order to 
understand the reasons for such tendency, we took ad- 
vantage of the LEIP and calculated partial contributions 
to the total exchange constants. 

We found that as in NaTiSi2 02, in V-pyroxenes the 
most important contribution is the direct exchange in- 
teraction between t2g orbitals. For LiVSi202 Jt2g-t2g — 
59.6 K. The t2g — eg partial exchange is FM, since it oc- 
curs from occupied to empty orbital. Numerically this 
contribution is defined by Eq. and in our LSDA-I-U 
calculation for LiVSi202 Jt2g-eg = —9.9 K. The rest, 
0.6 K (total exchange for LiVSi2 02 is 50.3 K) is the par- 
tial exchange between eg orbitals, which should be com- 
pletely empty in the purely ionic picture, but which still 
are partially occupied in the real band structure calcula- 
tion because of the hybridization with oxygens. 

Going from Li to Na, the t2g — t2g partial exchange 
weakens because of the increase of the distances and local 
distortions of the octahedra. For NaVSi202 Jt2g-t2g — 
27.5 K, Jt2g-eg = -13.3 K and J^^-e, = 0.8 K. The 
direct calculation of hopping integrals corroborates this 
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FIG. 8: (color online). Total and projected DOS for 
NaCrGe206 obtained in LDA+U calculation. Fermi level is 
zero energy. 



interpretation. Thus e.g. for LiVSi206 and NaVSi206 
the ratio of t2g — t2g exchange integrals J^^* /Jt^"" = 2.17, 
which is close to the squared ratio between effective t2g — 
t2g hopping integrals (i^Vt^")^ = 2.15. 

The total values of intrachain interaction obtained 
within LMTO method J = 50.3 K (LiVSizOg) and 
J = 15.0 K (NaVSi206) agree with those calculated 
in pseudo-potential code from the total energies: J = 
51.3 K (LiVSiaOg) and J = 13.8 K (NaVSiaOg). In 
addition it matches the exchange constants obtained in 
Ref. IIO, where the fitting of experimental curves of mag- 
netic susceptibility to theoretical one, obtained using 
Quantum Monte Carlo simulations on a cubic lattice, was 
performed. Note however, that Pedrini et aZ.jl^ used dif- 
ferent definition of Heisenberg Hamiltonian, and their ex- 
changes must be divided by factor two to compare with 
ours. In addition it should be mentioned that the real 
geometry of pyroxenes includes frustrating diagonal ex- 
change paths depicted in Fig.[2la) as dashed Hues. Those 
contributions were not taken into account in Ref. l40l 



VI. Cr-BASED PYROXENES {S) 



from the DOS plot obtained for NaCrSi206 in LSDA-hU 
calculation and presented in Fig. [5] Quite similar situ- 
ation occurs in other Cr-based transition metal oxides: 
one may observe significant mixing between O— 2p and 
Cr— 3d states at the top of the valence band [42|, H^. 
This feature of the electronic structure is important for 
the magnetism of Cr-based pyroxenes, since it increases 
super-exchange interaction via 0—2p states. 

Experimentally LiCrSi206 and NaCrSi206 were found 
to be AFM with small Neel temperature, 11 K and 3 K 
correspondingly; NaCrGe206 is FM with Tc = 6 K. 

In our calculations LiCrSi206 was indeed found to be 
AFM with J = 3.4 K. NaCrSi206 shows rather small 
intrachain exchange, less than 1 K. This number is ob- 
viously lying beyond the precision of the calculation 
scheme [4J]. One may just note that this compound is 
on the borderline between FM and AFM ordering. In 
agreement with the experiment, intrachain exchange in 
NaCrGe206 is FM and more pronounced, J = —5.2 K 
(—7.9 K in pseudo-potential calculation). 

More detailed investigation of the partial contributions 
to intrachain exchange shows that t2g — eg exchange in- 
teraction in all Cr-based pyroxenes is almost the same as 
in V-based ones: J|f_e^ -10.5 K. The AFM t2g - t2g 
contribution, in contrast, is significantly reduced in Cr 
pyroxenes. This seems to be the joint effect of the in- 
crease of the TM-TM distances and on-citc Coulomb re- 
pulsion U and the decrease of the charge-transfer energy 
A and radii of localization of d— electron rd, as it was 
discussed in the end of Sec. [Ill 

The largest t2g—t2g exchange was found in LiCrSi206: 
Jt2g~t2cj =13.4 K, which has shortest Cr — Cr distance 
(D = 3.06 A). For NaCrSi206 the t2g — t2g exchange is 
smaller, Jtag-ta, =7.7 K, which is presumably related to 
longer Cr~Cr distance {D = 3.086 A). Thus, according 
to our calculations the FM t2g — eg exchange almost com- 
pensates AFM t2g ~t2g contribution in NaCrSi206. This 
tendency continues for NaCrGe206, where because of 
larger Ge and Na ionic radii the Cr — Cr distance exceeds 
its maximum value {D = 3.142 A), and the FM t2g — Cg 
contribution becomes dominating, since -t^ = 3.1 K. 



Generally the electronic properties of Cr-based {d^ con- 
figuration, S = 3/2) pyroxenes are similar to those of 
other pyroxenes: they are also insulators, but with larger 
band gap ^3 eV. This is mainly because of the fact that 
spin-majority t2g sub-shell becomes fully occupied and 
goes lower in energy. More important is that the charac- 
ter of the valence band in Cr-based pyroxenes is different. 

Both the Ti- and V- based pyroxenes were found to be 
Mott-Hubbard insulators, with the gap separating occu- 
pied and unoccupied transition metal 3d bands. Going 
from the left to the right in the same row of the periodic 
table, the charge-transfer energy decreases, 0-2p states 
shift up [iHl and get to the same energy region as the 
lower transition metal Hubbard band. This is clearly seen 



A. Changes of the electronic structure going from 
Si to Ge on the example of Cr-based pyroxenes 

The Cr— based pyroxenes were the only pyroxenes with 
available and reliable crystal structure to study the effect 
of Si Ge interchange in these compounds. We studied 
only these, but argue that the general character of the 
changes in the electronic structure should be the same 
for other pyroxenes. 

There are two possible mechanism, which may deter- 
mine the difference in the electronic properties of Ge and 
Si based pyroxenes. First of all Ge— 4p states have larger 
spatial extensions than Si— 3p, this may result in larger 
covalcncy effects. Secondly the size of Gc^+ ions is much 
larger than Si^"'' (the ionic radius -Rg^4+ = 0.39 A, while 
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FIG. 9: (Color online). Cr-3d DOS for NaCrSi206 and 
NaCrGegOe obtained in LDA. It illustrates the decrease of 
the t2g — eg crystal-field splitting with substitution Si by Ge. 
Fermi level is zero energy. 



= 0.24 A). Therefore, different structural distor- 
tions are expected for Ge and Si pyroxenes. 



In order to understand how important is the effect of 4p 
or ip states, we performed the calculations, where in the 
NaCrGe206 Ge ions were substituted by Si. In this case 
all the distances in "artificial NaCrSi2 06" were taken the 
same as in the parent NaCrGe206. AH the radii of MT- 
spheres were also taken the same as for NaCrGe206- The 
results show minor difference in shape and positions of 
Cr— d bands compared to NaCrGe206. Therefore one 
may conclude that the origin of the difference in the elec- 
tronic properties between Si and Ge pyroxenes is mostly 
due to different local distortions. 

The average Ge~0 bond distance {dce-o ~ 1-74 A in 
NaCrGe206 ) is larger than Si — O one {dsi~o = 1-63 A 
in NaCrSi206 )• Since the TM-0 distance is mostly de- 
termined by the strong ionic TM-0 bond and hence by 
the TM and O ionic radii only, the substitution of Si by 
larger Ge leads to the bending of TM-O-TM bond. If 
in NaCrSi206 the angle Cr - O - Cr ^ 99.6°, then in 
NaCrGe206 this angle is 101.2°. This leads to an increase 
of the Cr—Bg bandwidth, while the Cr— t2g bandwidth is 
not changing, as one may see in Fig. [HI 

On the other hand, shorter Si — O bond distance re- 
sults in the larger shift of the 0—2p states to the lower 
energies, than in the case Ge. This should lead to an in- 
crease of the crystal-field splitting in TM— Sd shell going 
from Ge to Si. The direct center of gravity calculation 
shows that t2g — eg crystal- field splitting in NaCrSi2 06 
equals to 1.89 eV, while for NaCrGe206 it is 1.68 eV. 



We were able to find the crystal structure for only one 
Mn-based pyroxene: NaMnSi206 (mineral namansilite). 
Having one electron in eg sub-shell, Mn ion is Jahn- Teller 
(JT)-active. As a result the "distribution" of long and 
short Mn-0 bonds in MnOg octahedron is different from 
other pyroxenes. 

It is known that because of the anharmonic effects 
the JT distortion around Mn^+ ions tymcally leads to 



a local elongation of MnOg octahedra [45[. These elon- 
gated octahedra may in principle be packed differently, 
leading both to an antiferrodistortive (for example as in 
LaMnOa) and to a ferrodistortive (e.g. Mn304) ordering. 

In the edge-sharing geometry the largest gain in the 
elastic energy may be achieved, if the longest Me-0 bonds 
in neighboring octahedra are parallel to each other. The 
explicit expressions for the elastic energies of different 
distortion modes are summarized in Tab. 1 of Ref. |46| . 
Among various possible types of parallel JT distortions 
in the edge-sharing octahedra, the consecutively "twist- 
ing" geometry of pyroxenes selects the only variant: the 
distortion occurs along the 0-Mn-O bond connecting two 
neighboring octahedra (local z— axis). Only in this case 
the systems gains the energy for both neighbors. The 
local crystal-field due to the presence of "side-oxygens" 
forming short Mn-0 bond, directed only to the Si04 
tetrahedra, stabilizes single Cg electron on the Sz^ — 
orbital and therefore works in the same directions. 

This conclusion is in full agreement with the crystal- 
structure studies performed in Ref. IH, where the paral- 
lel 180° 0-Mn-O bonds along the local z— axis connect- 
ing neighboring MnOg octahedra were found to be the 
longest. 

In the LSDA-I-U calculations we indeed obtained that 
on all of Mn'^^ ions the single eg electron localizes on the 
same — orbital, directed along this longest Mn-0 
bond. Resulting magnetic moment equals to 3.69 the 
bang gap ~ 1.76 eV. The intrachain exchange calculated 
in LEIP procedure is AFM and ~ 3.4 K (from the total 
energy difference in pseudo-potential code we obtained 
that J = -0.1 K) 

The t2g — t2g and t2g — eg contributions to the total 
exchange are AFM and equals to 2.1 K and 6.1 K respec- 
tively. The situation with half-filled t2g — <2g sub-shells 
is the same as in the case of Cr, but t2g — eg exchange 
will now predominantly go from the occupied t2g to occu- 
pied Cg orbitals, as shown in Fig.[3Jb). This contribution 
should be AFM, as follows from Eq.(I3]). The Cg - Cg 
contribution is about —4.8 K FM. 



VIII. Fe-BASED PYROXENES {(f) 

The Fe'^''" ions have the simple configuration d^{t2ge^), 
without any orbital degeneracy, thus the exchange is uni- 
form along the chain. Using LEIP formalism we obtained 
that the intrachain exchange is AFM both in LiFeSi206 



(J=15.8 K) and NaFeSisOg (J=15.9 K). This is con- 
sisted with neutron measurement performed in Ref. for 
LiFeSi206, but disagrees with the experimental findings 
for NaFeSi2 06 presented in Ref. [l^, where Fe-ions were 
claimed to be ferromagnetically ordered in the chain. 

The t2g — t2g partial contribution to the exchange in 
Fe pyroxenes was found to be AFM and of order of that 
for Cr-based ones: 0.5 K and 1.4 K for LiFeSi206 and 
NaFeSi206 respectively. The t2g — Cg exchange is also 
AFM since it goes from the occupied t2g to occupied Cg 
orbitals, as in the case of NaMnSi206. It equals to 9.4 K 
and 7.0 K, respectively. 

The situation with Cg — Cg contribution here is a bit dif- 
ferent from the case of NaMnSi2 06. The TM-O-TM an- 
gle equals to ~ 97.1° in NaMnSi206, while in NaFeSiaOg 
and LiFeSi206 it is larger: ~ 100.5° and - 99.0° re- 
spectively. This may be significant, since in deviating 
from 90°, small FM contribution given by Eq. ^ is 
quickly outbalanced by much stronger AFM exchange 
of order /Udd- Detailed analysis shows that the bor- 
derline, where the mutual compensation occurs, is about 
97° [43]. Thus, one may expect that in Fe-based pyrox- 
enes Cg — Cg exchange will be already AFM. Note however 
that such a strong angle dependence is not a common fea- 
ture of all the 90° contributions described in Sec. [Ill For 
instance, the important for the cromatcs t2g~eg exchange 
is almost angle independent. 

Our LEIP calculation confirms this qualitative discus- 
sion. The Cg — Cg contribution in Fe pyroxenes is indeed 
AFM and equals to 5.9 K and 7.4 K for LiFcSi206 and 
NaFeSi206 respectively. 

Thus, from a theoretical point of view it is hard to 
expect that the intrachain exchange would be FM, as 
the authors of Ref. [H proposed. We see that the ex- 
change constants and their partial contributions obtained 
in LEIP formalism agree with the theoretical conclusions 
based on the perturbation theory. However, they are 
obviously too large to describe experimentally observed 
paramagnetic Curie temperatures: O = —25.8 K for 
LiFeSi206 H and 6 = -46 K for NaFeSisOg 

In contrast to other pyroxenes the results of total en- 
ergy calculations in the case of Fe-based pyroxenes differ 
from those obtained in LEIP formalism. The solution 
with AFM ordering in the chain was found to have low- 
est total energy for both Fe-based pyroxenes, and the 
intrachain exchanges extracted from comparison of AFM 
and FM total energies are 7 K and 8.5 K for LiFcSi2 06 
and NaFeSi206 respectively. The difference with the 
LEIP results appears because of the strong polarization 
of oxygen 2p-shell. This contribution is not explicitly 
(only via hybridization effects) taken into account in the 
present LEIP formulation, since both the Green func- 
tions used in Eq. dH) and defined in Eq. ([9]), and the 
potential correction Aij are calculated only for transi- 
tion metal 3d— states. The values of exchange constants 
obtained from the total energies agree much better with 
the experimental values of the Curie temperature. 

The detailed explanation of the origin of the diffcr- 
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spin-polarization on oxygen ions in Fe-based pyroxenes. 

ence in the results for the exchange constants obtained by 
LEIP and from the total energies for Fe pyroxenes, is the 
following. According to the crystal structure of pyrox- 
enes all the oxygens surrounding every transition metal 
ion can be divided on two categories. There are two "side- 
oxygens" finking TMOg octahedra and Si(Ge)04 tetrahe- 
dra, and four oxygens connecting neighboring octahedra, 
as shown in Fig. El^b). In the case of AFM intrachain 
coupling, those four oxygens are situated in the compen- 
sated field and will not have any polarization. Two other 
oxygens may be magnetized if the exchange field is large 
enough. For the FM chains both types of oxygens will 
be polarized. Such a polarization will lead to an addition 
magnetic contribution in the total energy calculation. 



This effect is especially important and pronounced for 
the pyroxenes, based on the late transition metal ions. 
Those ions have large magnetic moment and hence could 
easier magnetize oxygens. Even more important is the 
other factor: the non-zero magnetization of oxygen im- 
plies strong mixing of magnetic d— states of TM with one 
of 0—2p orbitals of the same spin, and therefore it will 
be much stronger for TM ions with partially filled Cg sub- 
shell. For Fe-based pyroxenes with half-filled Cg shell this 
effect is rather large: the average magnetic moment on 
oxygens, lying between TMs, is 0.17 fiB for LiFeSi206 
and 0.16 fiB for NaFcSi206. The oxygen moments are 
directed parallel to those of Fe. It is clearly seen from 
Fig. [ini only hoppings with antiparallel spin are possible 
between O— 2p to Fe— Sd shells. 

Additional polarization of 0—2p shell in FM solution 
will lower its total energy. Qualitatively one may consider 
it as following. In the AFM case both spin up and down 
electrons would hop from oxygen to different iron ions, 
while only spin down electrons will be active ones in FM 
situation, as shown in Fig. [TOl Because of the Hund's 
rule coupling on oxygen, the latter (FM solution) will 
gain more energy. 

It is easy to estimate the lowering of the FM total 
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energy because of the oxygen polarization recalling that 
the gain in magnetic energy in the band theory equals to 

Eo - -/MV4, (10) 

where / is the Stoner parameter defined as the second 
derivative of the exchange-correlation energy with re- 
spect to the magnetic moment M: 

I - ~2d''E,Jd''M. (11) 

The Stoner exchange constant / for the oxygen 2p-shell 
was estimated to be / = 1.6 eV in Ref. |4^, using 
constrained LDA calculations. 

Recalculating the total energy contribution given by 
Eq. pU]) in terms of the exchange constants for Heisen- 
berg Hamiltonian and taking into account that there are 
four polarized oxygens per two formula (without two 
"side-oxygens") we obtained that the exchange correc- 
tion due to the oxygen polarization is 10.7 K and 9.5 K 
for LiFeSi206 and NaFeSi206 respectively. The differ- 
ence between these values and those number obtained in 
LEIP formalism is a little bit larger than the exchange 
extracted from the total energy calculations. This seems 
to be because we did not account for the changes in the 
one-electron energies (relaxations), which tends to com- 
pensate the magnetic corrections given in Eq. (fTO|) . 

Finally, it is important to note that the contribution, 
which comes from the polarization effects on oxygen, is 
not a unique feature of pyroxenes. It was recently found 
to be crucial for the description of the magnetic inter- 
actions in Li2Cu04 [i^. Authors of Ref. originally 
interpreted this effect in the manner of initial ideas of 
Heisenberg as the exchange between Wannier functions 
(constructed out of TM-d and 0-2p orbitals) due to the 
Coulomb interaction. It is interesting to note that the 
final expression for the correction presented in Eq.(22) of 
Ref. \^ is exactly the same as ours written in Eq. (fTO|) . if 
one uses substitution: (3 = M{0)/2. 

IX. INTERCHAIN EXCHANGE 

In addition to the intrachain exchange interaction we 
calculated for Fe-based pyroxenes the interchain ones. It 
is especially important for the analysis of apossible origin 
of multiferroicity in i<"e-based pyroxenes [l|. The results 
reflect the general tendency for interchain exchanges in 
the whole pyroxenes family. There are two types of inter- 
chain paths, as depicted in Fig. [21^ a): some of TM ions 
in different chains are connected via two Si04 (GCO4) 
tetrahedra (J2), while the others - via only one (Ji). 

First of all. we found that the interchain exchange is 
much weaker than the intrachain one. Then, we found 
that for both LiFeSi206 and NaFeSi206 J2 is, indeed, 
approximately two times larger than Ji: jf ' = 1.9 K, 
J^' = 3.4 K, » = 0.8 K, J2^° = 1.6 K. In addition 
one may see that the change of Li ions to Na leads to 
a significant decrease of interchain exchange interaction. 



It is connected with larger ionic radii of Na (in 6-fold 
coordination i?j^j+=0.76 A, while i?^^+=1.02 A). 



X. CONCLUSIONS 

In this paper, using ab-initio LSDA+U calculations, 
we carried out the detailed analysis of electronic struc- 
ture and especially of the exchange interaction in a broad 
class of quasi-one-dimensional magnetic systems - pyrox- 
enes containing transition metal ions. We analyzed the 
systematics of the exchange interactions and compared 
it with the simple rules following from the perturbation 
theory. 

The exchange interaction due to direct hopping be- 
tween t2g orbitals (which we call "direct exchange inter- 
action" ) dominates in early transition (Ti,V) metal-based 
pyroxenes. In NaTiSi206 single d— electron localizes on 
the zy— orbital, which results in a strong AFM coupling 
in one out of two Ti-Ti pairs. This leads to the dimeriza- 
tion of TiOg chains and the formation of a spin gap on 
singlet Ti-Ti dimers. 

The nonuniformity of the exchange interaction disap- 
pears in (P pyroxenes, such as (Li,Na)VSi206, where elec- 
trons occupy zx and zy orbitals. This type of orbital 
filling is determined by a common structural feature of 
all pyroxenes (except Jahn- Teller ones): the presence of 
two short TM-0 bonds perpendicular to a local z— axis 
(directed along single 180° O — TM — bond connecting 
two neighboring TMOe octahedra). 

With an increase of the number of d— electrons AFM 
t2g — t2g direct exchange interaction is gradually sup- 
pressed due to an increase of the on-cite Coulomb inter- 
action ([/), and the t2g — t2gl&g super-exchange contri- 
butions (most of which are FM) starts to play a more 
important role, especially because of the decrease of the 
charge-transfer energy, which enters the corresponding 
expressions for the exchange in Eqs. in denomi- 

nator. In Cr-based pyroxenes (d^) the AFM t2g - t2g 
exchange interaction is nearly compensated by the FM 
t2g — exchange. The crystal structure makes a fine- 
tuning of these contributions, and NaCrGe206 with 
largest Cr — Cr distance turns out to be FM. 

In NaMnSi206, Mn^+ is a Jahn- Teller ion. We show 
that the consecutive "twisting" geometry of pyroxenes 
and the tendency to gain an elastic energy result in the 
ferro type of distortion with the local z-axis becoming 
the longest one. 

Iron-based pyroxenes were found to be AFM, with the 
large oxygen polarization (/iq ~ 0.2/1^). This effect is 
usually ignored in the literature [l3| . However, we show 
that it considerably reduces (by ~ 2/3) the AFM super- 
exchange interaction. The total energy of the FM state 
is lowered, because of the gain in Hund's rule coupling in 
O— 2p shell, which turns out be polarized in this state. 

The interchain exchange coupling was estimated for 
Fc-bascd pyroxenes and found to be small and AFM. 
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The presented results give a rather complete descrip- 
tion of exchange interactions in magnetic pyroxenes, and 
it may serve as a basis for the further study of magnetism 
and other properties, including multiferroicity, of this in- 
teresting class of materials. 
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